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O : ^^^^'^^ developn^ents in the physics of low density trapped gases make it worthwhile to verify old, 

' well known results that, while plausible, were based on perturbation theory and assumptions about 



(N 



(N 



X 



pseudopotentials. We use and extend recently developed techniques to give a rigorous derivation 



^ ' of the asymptotic formula for the ground state energy of a dilute gas of fermions interacting 

■ with a short-range, positive potential of scattering length a. For spin 1/2 fermions, this is £ ~ 

■ + {h^ /2m)2-KN ga, where ij" is the energy of the non-interacting system and g is the density. A 
similar formula holds in 2D, with ga replaced by g/\ \u{ga^)\. Obviously this 2D energy is not the 
expectation value of a density-independent pseudopotential. 



> 

o 
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Q ; I. INTRODUCTION 

(N 

Tij- \ The leading asymptotics for the ground state energy of a dilute gas of fermions, interacting with a positive, short 
, range pair potential, was derived years ago by several approximate methods . Indeed, the leading correction 

beyond the ideal gas formula is no different for fermions than for bosons, except for the fact that all bosons interact 
■ with each other whereas the spin-up fermions effectively interact only with the spin-down fermions and not with each 
' other. In this sense, the formula goes back to Lenz Q who derived the energy formula for bosons by assuming that 
each particle interacts with N — \ fixed particles that are well spaced from each other. 

Thus, we expect that the ground state energy per unit volume, e(gi|, qi), for spin-up particles and spin-down 
particles of mass m in a box of volume V (in the usual thermodynamic limit in which V ~* oo and = N-^/V and 
Ql = Ni/V are fixed) is, asymptotically, 

e{Q], Qi) = ^^(^'^^)^^^ (^'f ^ + ^) + a + higher order in [q^^qi) , (1) 

where a is the two-body (s-wave) scattering length of the pair potential v. Under the assumption that the total 
density g = g-^ + g i la fixed, this formula indicates that at low density the absolute ground state has spin zero, i.e., 
= = f/^- 

The corresponding low density formula for bosons contains only one kind of density and is 

e(g) = Aiiag^ + higher order in g . (2) 

This asymptotic formula was proved rigorously in Q . 

It is customary, nowadays, to regard (^) as coming from a pseudopotential ^87rai5(a;i — Xj), and that is certainly 
a useful shortcut to obtaining current results. But this has to be justified mathematically, and that is the purpose of 
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this paper. Several issues of physical interest are involved, which make it not totally obvious that the pseudopotential 
approach is beyond need of justification. 

• The availability of a good variational function 5* is important in theoretical physics, one that correctly displays 
the correlations of physical interest and whose energy (which is necessarily an upper bound to the 
ground state energy) can be computed without resort to uncontrolled approximations. It should also give the 
correct energy to the desired accuracy. A good example is the BCS function of superconductivity. In the boson 
problem one would think of a Bijl-Dingle-Jastrow function J = Y[i j gi^i ~ Xj) but it has not been possible, as 
far as we know, to carry out the energy calculation without making assumptions. The correlations are subtle 
(even if they are physically clear) emd have to be treated carefully, and the required rigorous bosonic upper 
bound was finally found by Dyson 6] but by using a non-bosonic variational function. In the fermionic case 
considered here we use a function of the form = S ■ J, where S" is a Slater determinant. While this ^E* looks 
simple, the calculation of an upper bound of the required accuracy ^ occupies half of this paper, and one 
cannot say that this is a simple calculation. 

• The source of difficulties in the boson problem is the subtlety of the correlations, which constitute the entire 
energy. In the fermionic case, on the other hand, we are looking for a tiny correction to a dominant free- 
particle kinetic energy, but this contribution, especially for hard-core potentials, does not come from a small 
perturbation. It is not obvious that hard core collisions do not create energy changes by perturbing the Fermi 
surface. 



The pseudopotential idea, while attractive, does have the drawback that it cannot be right in two dimensions 



(2D). The quantity ag'^ for 3D bosons is replaced by /\^n{ga'^)\, as predicted by [2|,I3 and proved in (Note: 



The scattering length can be defined in 2D as well as 3D. See |9j.) Consequently, the pseudopotential will have 
to depend on g. Thus, for fermions the energy to leading order in ga^ ought to be 

e{Ql,Ql) = ^27r (^?^ + g^) + ^ | ij^(g^2)| ^Tgi + higher order in {g^, gi) . (3) 

We will prove formula Q as well. Fermions in two dimensional layers are physically interesting and it can be 
useful to have Q proved rigorously. 

• In addition to the pseudopotential approach there is also the approach of summing diagrams Q|, which leads 
to many terms beyond the two in Nevertheless, it has to be admitted that expansions, especially where 
hard-core potentials are concerned, may have convergence or other difficulties. Effective field theory methods 
have also been used successfully but with similar concerns. Therefore, rigorous confirmation is much to be 
desired and we provide it here. 

In the following we do everything in 3D until Section IVII where we explain the modifications necessary for 2D, 
some of which are not trivial. This is done in order to make the essential ideas as clear as possible. In a forthcoming 

n 

paper 12], the natural generalization of to positive temperature states will be proved. 
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II. MODEL AND MAIN RESULTS 

In units in which f? jlm — 1 (which wiU be used henceforth), and with A = V^, the Hamiltonian is given by 

N 

H = Y,-^^+ vi^^-^3)> (4) 

i=l l<i<j<N 

acting on anti-symmetric functions of N space-spin variables, i.e., functions in /\^ i^(A; C^). Here q> 1 denotes the 
number of spin states. The particles are confined to a bounded region A, which we choose to be a cube of side length 
L and volume V = (or in 2D). We choose Dirichlet boundary conditions for the Laplacian, i.e., '5 = when 
any Xi is on the boundary of the cube. 

Since H is spin- independent, we can specify the number of particles of each spin, iVi, N2, ■ ■ ■ , Nq with N — Nj. 
The wave function ^ is then a function of the N coordinates xi, . . . , xtv, without mention of spin at all, but with the 
requirement on 'I' that it be antisymmetric separately in the first A^i variables, the second N2 variables, etc. To avoid 
needless notation we will give our proofs for q = 2 but will state the main Theorems ^ a-nd El for general q. 

The pair potential v{x) is assumed to be positive, radial, and of finite range Rq. It then has a finite and positive 
scattering length a. The scattering length can be defined as follows: if fix) is the unique solution (see 9] for a full 
discussion) of the zero-energy scattering equation 

-A(p(x) + ^v{x)ip{x) = (5) 

subject to the boundary condition lim|2,|^oo fix) = 1, then a is given by a = lim|^|^oo |2;|(1 — fix)). Note that we do 
not assume v to be integrable; our results also apply to the case of a hard core. Note also that for a pure hard-core 
interaction, the scattering length is equal to the range. 

There is no need (apart from simplicity) to assume that the potentials w(a;) between different groups of particles 
are the same, thereby allowing the Hamiltonian to be 'spin-dependent'. Thus, we could take the pair potential to 
be Vij(x), (with 1 < i,j < q) between groups i and j, with corresponding scattering lengths a^j. In this way, the 
quantity agiQj in Theorem ^ would be replaced by aijgiQj. Our proof would still go through with obvious trivial 
changes. 

Our main result concerns the ground state energy EQ{{Ni}, L) of H, in the thermodynamic limit i — > 00 with 
Qi = Ni/L^ fixed. It is well known that for systems with short range interactions the limit of the energy density, 
Eoi{Ni}, L)/L-i exists and is independent of boundary conditions |l3lll4|. 

Theorem 1. Fix Qi = Ni/L^ for 1 < i < q and g = Qi, and let EgHNi}, L) denote the ground state energy of H 
with the appropriate antisymmetry in each of the Ni coordinate variables. Then, for small g, 

hm ij£;o({A^J,i) = ^(67r2)'/=^^^,p + 87ra ^ g, g, + ag'eig), (6) 

with —const, iag^^^)^^^"^ < e{g) < -l-const. iag^^^)^^^ . 

The constants in the bounds on eig) depend on the interaction potential only through the dimensionless ratio 
Ro/a. We could, in principle, display the explicit dependence on Ro/a. By cutting off an infinite range potential in an 
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appropriate way, this would allow us to extend the result (with different bounds on s{g)) to infinite-range potentials 
with finite scattering length. 

The analogous theorem in 2D is the following. 

Theorem 2. Fix Qi — NijL? for I < i < q and g = Qi, and let EQ{{Ni}, L) denote the ground state energy of H 
with the appropriate antisymmetry in each of the Ni coordinate variables. Then, for small g, 

with -const. I ln(a2g)|"i/i'' <£(£») < +const. | ln(a2gi)|-i ln|ln(a2e)|. 

For simplicity we consider only q = 2 henceforth, i.e., the spin i case. The extension to general g > 2 is straightfor- 
ward. We introduce the following convenient notation. For A^i + A^2 = N, let X — {xi, . . . , xni) and F = (j/i, . . . , j/at^) 
stand for the collection of spin-up and spin-down particle coordinates, respectively. The Hamiltonian can then be 
written as 

H = -Ax - Ay +VXX +VYY +VXY, (8) 

with Ax = = J2^'=l Ay = = J2f=i ^li^ ^xx = E,;<j ^i^'i ~ and vxv = E^j - Vj)- It acts on 
the Hilbert space of square-integrable functions that are antisymmetric in the X and in the Y variables. 



III. OUTLINE OF PROOF 



Before presenting the proof of Theorems ^ and [21 in full detail, we give a short outline. We first concentrate on the 
three-dimensional case. The proof is split into two parts, the upper and lower bounds to the ground state energy. The 
upper bound, given in Section llVI uses the variational principle. The idea is to construct a trial wave function that 
shows the features one would expect the true ground state to have, but is at the same time sufficiently simple to make 
it possible to compute a good upper bound on the expectation value of the Hamiltonian. For this latter purpose we 
find it necessary to choose a function that confines the particles to small boxes, separated from each other to avoid 
interaction between different boxes. These boxes must not be chosen too small, however, to ensure that the finite 
size effects are negligible compared to the leading term in the interaction energy. Since this latter energy is rather 
small, we are forced to have a large number of particles in each of the boxes. This makes it impossible to control the 
norm of our trial wave function, and hence we must carefully take into account cancellations between the ratio of the 
expectation value of the Hamiltonian and the norm of the trial wave function. 

For the lower bound to the energy, given in Sectional the first essential step is to replace the 'hard' interaction 
potential v[x) by a 'soft' one, W(x), at the expense of using up some positive kinetic energy. This idea goes back 
to Dyson 6], who computed a lower bound for the ground state energy of a hard-sphere Bose gas. Only the high- 
momentum part of the kinetic energy is dispensable, however, since the low-momentum part is needed to fill the 
Fermi sea. In Lemma 0] below we prove such a bound, using only momenta bigger than a certain cutoff. With the 
soft potential W(x) one can then hope to proceed with some sort of rigorous perturbation theory to obtain a lower 
bound to the energy. Indeed, we prove two a priori bounds, one on the one-particle density matrix of the ground 
state, showing that it is close to the projection onto the Fermi sea, and another one on the number of particles whose 
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distance to their nearest neighbor is small. These bounds can be used to show that the ground state expectation of 
W{x) has the anticipated value. 

The necessary modifications of our proofs for the corresponding result in two dimensions, Theorem |21 are sketched 
in Section IVll 

IV. UPPER BOUND TO THE GROUND STATE ENERGY 

We start by collecting some properties of the solution to the 3D zero-energy scattering equation The proofs 
can be found in the appendix of P|. The solution to (jSJ, 1^5(2;), is a radial function and satisfies 



< (p{x) < 1, and hence a > 0. 



• ip{x) is subharmonic on (i-e., Aip[x) > 0, see [isl'). Aip{x) is a positive measure which is zero for |x| > Rq, 
and Jjj3 Aip{x) (Px = Aira. 

• '^{^) ^ 1 ^ and ip{x) = 1 — a/\x\ for |a:| > Rq. 

• !\x\<R + \vMx)?) d\ = 47ra(l - a/R) for R > i?o- 

These properties will be useful both for the upper bound given in this section and the lower bound given in the 
next. 

For the upper bound, it will be convenient to localize the particles into small boxes with Dirichlet boundary 
conditions. The number of particles in each box will be large for small p, but finite and independent of the size of the 
large container V. Let the side length the small boxes be £. If we place these small boxes a distance Rq from each 
other, then there will be no interaction between particles in different boxes. We then want to put n = gi{£ + Rq)^ 
spin-up particles into each box, and likewise m — g2{£ + Ro)^ spin-down particles. Since gi{£ + Rq)^ need not be an 
integer, however, we will choose 

n = gi{i + Rof + ei and m = g2ie + Rof + e2, (9) 

with < £i,£2 < 1 chosen such that n and m are integers. We then really have too many particles, but this is 
legitimate for an upper bound, since the energy is certainly increasing with particle number. We thus have 

lim l^EoiNi,N2,L) < ^ Eo{n,m,e), (10) 

with n and m given as in ©. This bound holds for all choices of the box size £. 

We will now derive an upper bound on the ground state energy of n spin-up and m spin-down particles in a cubic 
box of side length £, for general n, m and £. We take as a trial state the function 

^{X,Y) = DniX)D^iY)GniX)GmiY)F{X,Y), (11) 

where Dn{X) denotes the Slater determinant of the first n eigenfunctions of the Laplacian in a cubic box of side 
length £, with Dirichlet boundary conditions. (In the case of degeneracy, any choice will do.) Moreover, 

Gn{X)= n 5(2^.-2:,), (12) 

l<«<j <n 
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with < g{x) < 1, having the property that g{x) = for |a;| < s and g{x) = 1 for \x\ > 2s, for some s > 2Ro to be 
chosen later. We can assume that |Vfif(x)| < const. for some constant independent of s. Finally, 

n rn 

F{X,Y)=Y[Y[f{x,-yj), (13) 
j=i j=i 

with f{x) = <fi{x) /{1 — a/R) for \x\ < R and 1 otherwise. Here (fi{x) denotes the solution to the zero-energy scattering 
equation, and we assume that R > Rq, which guarantees that / is a continuous function. Moreover, we assume 
2i? < s. By the variational principle, 

Eo{n,m,i)<^±^. (14) 

Since ^ vanishes whenever two particles of the same kind are closer together than the range of the interaction, we 
have 

(*|ir|*) = (*| - Axl*) + (*| - Arl*) + 

In evaluating the kinetic energy, we use partial integration and the fact that Z)„(X) is an eigenfunction of — A^. Let 
the corresponding eigenvalue (namely the sum of the lowest n eigenvalues of the Dirichlet Laplacian) be denoted by 
E^{n,£). Then 

(*|-Ax|*) = S°(n,^)(*|*) 

+ j Dr,{X)'\VxGn{X)F{X,Y)\^D,^{YfGm{YfdXdY. 

Here we denoted dX = HlLi ^^'^ ~ Hjli ^'^^ short. In the second term, we use the Schwarz inequality 
to deduce (for some e > to be chosen later) 

\VxGn{X)F{X,Y)\^ < {l+e)\WxFiX,Y)\'G^iXf 

+ (! + £-!) F{X,Yf\VxGn{X)\^. 

Proceeding in the same way for the kinetic energy of the F-particles, we thus get the upper bound 

{^\H\^) < I + (1 + £)ll + (1 + £-^) III, (15) 

with 

1= [E°(n,^)+ i;D(m, £)](*!*), (16) 



H = 1 [\VxF{X, y)|2 + \VyF{X, Y)f + vxyF{X, Yf] Dr.{XfDmiYfGn{XfGm{Y)^ dX dY, (17) 

and 

HI = j [\VxGn{X)\^Gm{Yf + \VYGm{Y)fGn{Xf] F{X, YfDr,{XfDm{Yf dX dY. (18) 
The positivity of vxy has been used here. We shall now bound these three terms, when divided by (^'I^E'), separately. 
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We start with I. We may consider the sum of the n lowest Dirichlet eigenvalues as a Riemann sum for the integral 



IpIS^f. 

PI ■P2 .P3 >0 

where we denote the Fermi momentum by fcp — {6n^ny/^ /£. It is then easy to see that 

E^{n,i) < ^(67r2)2/3!^ (l + const, n-^/^) , (19) 

and the exponent in the error term is actually optimal. We note that this bound shows that we must not choose £ too 
small in order to have an error term that is negligible compared with ag; more precisely, we need n ~ gi£^ ^ (a^p)"^. 
This will be fulfilled with our choice of £ below. 

Next we derive an upper bound on II. We are going to need the following lemma. 

Lemma 1. Let D„ {X) denote a Slater determinant of n linearly independent functions (j)a{x). For a given function 
h{x) of one variable, let ^{X) he the function = Dn{X) n"=i h{xi), and let M denote the n x n matrix 

M^p = [ 4,l{x)cf>p{x)\h{x)\'' d^x. (20) 



Then 

(i) The norm of ^ is given by ('I'I'I') = det Af. 
(a) For 1 < k < n, the k-particle densities of ^ are given by 



kj ($1 



/ \^{X)\^d^Xk+i---d^Xn = ■^[||/i(a;,)P(xiA---Axfe|M-i(g)---(g)M~i|a;i A---Axfc), 
' ' ' i=i 

where \x) denotes the n-dimensional vector with components (j)a{x), 1 < oi < n, and |a;i A • • • A x^) stands for 
the Slater determinant {k\)~^/'^ Eo-(~-'-)'^I^ct(i)) ® ■ ■ ■ ® |a^(T(fc)), cr denoting permutations. 

(Hi) If ^'^(X) — ^{X)k{xi) / h{xi) for some function k{x), then 

n 

^($^1$:) ^ (detA/) (Tr[i^Af-i]), (21) 
1=1 

where Tr [ • ] denotes the trace, and K is the n x n matrix 

Kc.p= [ rAx)Mx)\k{x)\^ d^'x. (22) 



The proof of this lemma is a straightforward exercise that we leave to the reader. Note that without loss of generality 
one can set h{x) ~ 1 in proving the lemma. Item (iii) is an immediate consequence of items (i) and (ii), noting that 
the left side of (|21|l . when divided by the norm of ^{X), is the integral of the one-particle density of ^{X) multiplied 
by \k{x)/h{x)\'^; we state it as a separate item for later use. 

Using Gn{X) < 1, we infer from this lemma that, for any fixed Y, 

J Gn{XfDn{Xf [\^xF{X,Y)\'' + \vxy\F{X,Y)\''] dX 
< I Dn{X)mVxF{X,Y)\^ + ^vxy\F{X,Y)\^] dX 

^TriKyMy'] j Dn{Xf\F{X,Y)\'dX. (23) 
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The n X n matrices Ky and My are given by (|20|l and (|22() . with 4>a{x) being the lowest n Dirichlet eigenfunctions 
of —A, and with h(x) — Y[j fix — yj) and 

\Hx)\^ = |v.n,/(^ - y,)\' + - y.)n,/G^ - Vi)^ 

respectively. Since Ky is a positive definite matrix, we have the bomid Tr KyMy^ < \\My^\\Tt Ky , where || • | 
denotes the matrix norm (i.e., the largest eigenvalue for hermitian matrices). To calculate Ti Ky, and to bound 
II My "'^11, we can assume that all the yj's are separated by at least a distance s, because the integrand of term II in 
(I17|l vanishes otherwise. Since s > 2R by assumption, we have in this case 

n 

\Hx)\^ = |v.n,/(2; - %)|' + hEjvi^ - y^m^fi^ - y,f = E " (^4) 

i=i 

with 

ax)^\Wf{x)\^ + \v{x)f{xf. (25) 
Hence, if Q^{x) denotes the one-particle density of Dn{X)^ we have 

n 

TlKy = Y.Sn*^iyj), (26) 

where * denotes convolution. 

r-l 



To bound ||My ||, we use the following: 
Lemma 2. Assume that \yi — yj\ > s > 2R for all i ^ j. Then 



1-My||< const. (^+n2/3£!K (27) 



Proof. Let q{x) = 1 — Y[j f{x — VjY — 0- Then, for any n-dimensional vector \b) with components b, 



{b\l-My\b)= j q{x)\Y,bMx) 



2 



Hence, the question about the largest eigenvalue of 1 — My translates into the question of how large the average 
potential energy for the potential q(x) can be for functions such as '^^ba^'aix) whose kinetic energy is bounded 
above by (const.) r?l'^t^'^ , i.e., the Fermi energy for n particles. 

Let Bj denote the ball of radius s/2 around j/j. Note that all these balls are non-overlapping by assumption. Also, 
since s > 2R, q{x) = if a; is outside all the balls. For a given function 4>{x), let (pj denote the average of (p{x) in the 
ball Bj. Moreover, let t]{x) = 0(a;) — (pj. By the Cauchy-Schwarz inequality (a + 6)^ < 2{a'^ + b^), we get the bound 

q{x)\(j){x)\^(fx<2 [ q{x)\ri{x)\^ d^x + 2\(j)j\^ [ q{x)d\. (28) 
Note that |^jf < &/{tts'^) /^^ \4>{x)\'^ d^x, again by the Cauchy-Schwarz inequality. Moreover, since s > R, 

q{x) d\ = / (1 - f{xf) d\ < (47r/3)ai?2. 



To obtain the last inequality, we used the definition of f{x) as well as the fact that (p{x) > max{l — a/|a;| , 0}, as 
explained in the beginning of this section. 
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Note that r]{x) is a function whose average over the ball Bj is zero. In other words, it is orthogonal to the constant 
function in Bj. Hence, using the fact that q{x) < 1 and Poincare's inequality 



qix)\'q{x)\^ d-^x < \r]{x)\'^ d-^x < const, \Vri{x)\^ d-^x. 

Bj JBj JBj 

In this last expression we can replace ri{x) by 4>{x), of course, since they only differ by a constant. Summing over all 
the balls Bj (and using that q{x) — outside the balls), we thus obtain that, for any function (t>{x), 



q{x)\4>{x)\ d X < const 



aR 



\(l){x)\^ d^x + / \V(l){x)\^d\ 



In the case in question, the kinetic energy of is bounded by const, ri^^^i ^. This finishes the proof of the 
lemma. □ 

Since < My < 1 as a matrix, this lemma implies that 

"^"^^ " " 1- 111 -My II - ^" " 1- const. [ai?Vs3 + n2/3(s/^)2] ' ^^^^ 

provided the denominator is positive. By inserting (|26|) and H29|) into l|23|) . we see that, for fixed Y with \yi — yj\ > s 
for all i ^ j, 

J Gn{XfD^{Xf [\^xF{X,Y)\'' + \vxYF{X,Yf] dX 

n „ 

<A^Y.^^* / D^{XfF{X, Yf dX. (30) 

To be able later to compare this expression (|30|) with (^l^*), we want to put Gn{XY back into the integrand. For 
this purpose we need the following lemma, which compares the integrals with and without the factor Gn{XY ■ 

Lemma 3. For any fixed Y, 

' D^{XfF{X,YfG^{XfdX 



> J Dn{XfF{X,YfdX (l-const.n'^^^\\M-^f{s/£)'^y (31) 
Proof. Since g{x) — 1 for \x\ > 2s, we have 

n 

G„(X)2>l-^0(2s-|a;,-a;,|). (32) 

i<j 

Here denotes the Heaviside step function, i.e., 9{t) = 1 for t > and 9{t) = for t < 0. To evaluate the integral 
of the second term in H32|l . we need the two-particle density of the state Dn{X)F{X,Y) for each fixed Y. By 
Lemma Q above, and the fact that f{x) < 1, this density, when appropriately normalized, is bounded from above 
by \\My'\\^g^'''^\x,x'), where gn'^^\x,x') denotes the two-particle density of the determinantal state Dn{X). In 
particular, by explicit computation one finds that this latter density satisfies the bound 

g^'^^Hx^x') < const. |x-a;'|2(n/£3)8/3 (33) 
for some constant independent of n and £. Hence we arrive at (|31f) . □ 
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We note that it is the n-dependence in inequality H31() that forces us to choose the particle number to be small and 
makes it necessary to localize the particles into small boxes. We also emphasize the importance of the exponent 5 
in (|31l) . which stems from the fact that the two-particle density vanishes as |x — x'p for x close to x' . Had we not 
taken this into account, we would get an error term of the order n?{s/t)'^ in (|31|l . Note that necessarily s > a, and 
hence this error would be huge if n 3> (a^ q)~^, which is demanded by H19(l . This also explains why it is not possible 
to treat F{X, Y) on the same footing as Gn{X) and GmiX), since the two-particle density only vanishes as \x — x'p 
for particles with equal spin, unlike the situation for particles of unequal spin. 

Let 

Bn= (l- const. n^/^Al{s/ef^ ^ , 
assuming that the term in parenthesis is positive. Applying Lemma |3 to inequality (|3()|l . we arrive at 
J Gn{XfD,,{Xf [|VxF(X,y)|2 + \vxYF{X,Yf] D„,{Y)G^{Y) dX dY 

n „ 

<A„B„^ / g'^*i{y,)D,^{YfDn{XfF{X,YfG^{YfGn{XfdXdY. (34) 
i=i 

Now we cannot bound * ^(y) independently of y by simply using the supremum of q^{x), since this number will 
be strictly bigger than n/^^, even in the thermodynamic limit. Instead, we repeat the above argument for the Y 
integration. We use |Gm(5^)| < 1, the F-analogues of Lemmanjand then LemmaOlto put GmiX)"^ back in. Here, it 
is important to note that now the XiS are separated by at least a distance s > 2R. In this way we obtain 

J GniXfD„{Xf [\VxF{X, + ^vxyF{X, Yf] Dra{Y)G^{Y) dX dY 

< A„BnBm J D^{YfD„{X)^F{X,YfGmiYfGn{XfTTKxM^UXdY. (35) 

The matrix Mx is the same as before, with Y replaced by X (and n replaced by m, of course), and Kx is the mx m 
matrix 



{KxU^ I MyrMy)l[fiy-^^?Qn*ay)d^. 



Using \f{x)\ < 1 and ||Af;^^|| < Am, which follows from Lemma |2 and the fact that the Xi's are separated at least by 
a distance s, we get the bound 

Tr KxM^' < A,„Tr Kx < An J g°(x)g° (2/)C(x - y) d\ d^y. (36) 

We recall the definition of f{x) and the properties of 'p{x) stated in the beginning of this section to calculate the 
integral j i{y)d^y = 47ra(l — a/R)^^. We then use this information to bound the last integral in H36|l. by using 
Young's inequality J^]. Thus 

Tr^xM^i < A„ (^j Q'^ixfd'^x^ ' [j (Pra{y?A^y^ ' 4™(1- a/i?)-i. (37) 
For the square of ^{x) we find 

3 

/ Q^X^f rf'^ = ^ E n (1 + ¥v.,,.) , (38) 

p,g a— 1 
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where pa denotes the components of the wave vector p, and the sums are over the n lowest eigenstates of the Dirichlet 
Laplacian, X^a=\ sin(paa;a)- From this expUcit expression it is easy to see that 

/2 
Q^S^f d\ < ^ (l + const, . (39) 

The same holds with n replaced by m. Eq. I|35|) thus implies the upper bound 

y G,,(XfD.^{Xf \^VxF{X,Y)\^ + \vxYF(X,Yf\ D,n{Y)G„{Y) dX dY 

< (^1^) ^^"""^ ^„A^B„g,„(l - a/R}-^ (^1 + const, n^i/^ + const. TO-^/^) . (40) 



The same bound holds, of course, with X and Y interchanged. We therefore have the upper bound 

II < (^|^)5!I^A„A„.B„B™(1 ~ a/R)-^ (1 + const, + const, . (41) 



£3 

It remains to bound the term III. Using |(7(a;)| < 1 we have that 



Xi ) I 



E E E |V.9(a:.-x,)||Vg(x,-a;fc)|. (42) 

^=1 i-J^i k,k^i,j 



Now, by Lemma n] the appropriately normalized fc-particle densities of Dn{X)F{X,Y) are bounded above by 
and Qn'''^^ satisfies 



My^W'' gn'''''\ where gn''"'^'' denotes the /c-particle density of DniX). In particular, Qn''"^^ is satisfies the bound 



e°'(3)(x,x',x") < const. (n/^3)3 

for some constant independent of n and £. Using the fact that \7g{x) is supported on the set < 2s, together with 
|V(7(a;)| < const. s~^, we obtain from Il42|l . for any fixed Y, 

J Dn{XfF{X,Yf\WxGr.{X)fdX 

2 r 

< const, {\\My^fn^/^s/£)^ + \\My^\fn{s/ef^ / D^iX)^ F{X,Yf dX. (43) 



Finally, to get a bound on III, we proceed as above, using H29|) (and the fact that the yj's are separated by a distance 
s) and Lemma O to put Gn{X)'^ back into the integral. Note, however, that it is enough to bound An and Bn 
by constants in this term. Assuming that n{s/£)^ is small, the second term in the parenthesis in (|43l) is negligible 
compared to the first term. The same bound applies to the case where X and Y are interchanged, and hence we 
obtain 

„3 

III < (*|*) const. (nS/3 ^ m^^^)—. (44) 

Collecting all the error terms obtained in Eqs. H19() . H41|l and H44() and inserting them into (|14|l and H15|l . we obtain 

3. 0.0/0 n5/3 + 771^/3 



Eo{n,m,i) < |(67r2)2/3!^^^(l + Cn-i/3 + Cm-i/3) 



nm 



(l + e + c\^ + {n + mf/^s/if + 1 + ^ + J- + {n + mf'\s/if 



Cs (n + m 
H 

e P 



l2 r 



(n + m)2/3(s/^)2l (45) 
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for some constant C > 0. In Ineq. (|45|l we have assumed smallness of all the error terms, i.e., that the terms in square 
brackets are small. This condition will be fulfilled, at low density, with our choice of R, s, n, m and £ below. 

The optimal choice of e in H45|) is given by — const, in + rnf' 1"^ j [(^ anrn) . Inserting this value for e we infer 
from 23 



+8Tra-^ {1 + C 



(l + C + (n + mf/^ is/if + ^ + ^ + ^ + in + ) 



g3/2^1/2 

+C(n + m)^/3 . (46) 

Eq. (|46|l is our final bound on the energy Eo{n, m, £). To apply this result in pU|) we have to insert the values @ for 
n and m. Recall that \n — gi{£ + Ro)^\ < 1 and |m — g2{£ + Ro)^\ < 1- We are then still free to choose R, s and £. 
We choose 

Inserting these values into H46|) we thus obtain, for small g, 

^Eo{n,m,£) < ^(6^2)2/3^^5/3^^5/3^ + 8TTagig2 + const. ag^ag^/^f^\ 

In combination with Eq. (|10|l . this finishes the proof of the upper bound. Note that the contribution to the error 
term that arises from the fact that Eo{n,m,£) has to be divided by {£ + 2Rq)^ and not £^ in (|l()|l . is of the order 
g^^^Ro/£ and, for this choice of £, is much smaller than ag^{ag^^^)^^^ when g is small. 

V. LOWER BOUND TO THE GROUND STATE ENERGY 
A. The Dyson Lemma 

We start with a generalization of a lemma of Dyson which bounds the hard potential v{x) (which may or may 
not contain a hard core) from below by a soft potential U{x), at the expense of using up some positive kinetic energy. 
In the following, f{k) denotes the Fourier transform of a function f{x), i.e., /(fc) = (27r)-3/2 / f{x) exp(ifc • x) (Px. 
Dyson's result was generalized in [5| and it is further generalized here by separating low from high momentum. In our 
application "low" will mean g^/'^ and "high" will mean 1/a. The analogous inequality for 2D is stated as Lemma [7| 
below. The proof of both lemmas is given in the appendix. 

Lemma 4. For R > Rq, let 9j^{x) denote the characteristic function of a hall of radius R centered at the origin, i.e., 
Or{x) = 1 if \x\ < R and — otherwise. Let x{p) ^2 a radial function, < x{p) ^ 1; such that h{x) = (1 — x){^) 
bounded and integrable. Let 

fnix) = sup \h{x -y)- h{x)l (47) 
\v\<R 



and 



wr{x) = ^fR{x) / fR{y)d\. (48) 
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Then, for any positive, radial function U(x), supported in the annulus Rq < \x\ < R, with Jjjg U{x) d^x = 47r, and for 
any e > 0, 

-Vx(p)0fl(a;)xb)V + ^vix) > (1 - e)aUix) - ^wnix). (49) 

Here, 9ji{x) is multiplication operator in x-space, whereas xip) multiplication operator in momentum space. Thus, 
'^xip)(^Rix)x{p)'^ fl*^ operator version ofV^, which is cut off in both configuration and in momentum spaces. 

The original Dyson lemma (as modified in 5]) has x{p) = 1 B^nd W]i{x) = 0, i.e., there is no cutoff. The cutoff xip) 
in H49I) essentially says that only the high momentum part of V is needed to give a good account of the scattering 
of two particles. The (relatively) low momentum part of V is not used in l|49() and is thereby saved for later use to 
give a good estimate of the part of the fermion kinetic energy needed to fill the Fermi sea. The price we pay for this 
luxury is the error term awji{x)/e, which does not appear in Dyson's lemma. 

Note that, by construction, either wb.{x) is bounded and integrable or else w^x) = oo for all x. If xip) = 1; 
then wj^{x) = 0, and hence we can set e = to recover the Dyson Lemma in '5], which says that for any 4>{x), 
^.|<«(l + iH^) " aU{x)Mxr) d\ > 0. 

Corollary 1. // j/i, . . . , yN denote N points in R"^, with \yi ~ yj] > 2_R for all i ^ j , then, as an operator on functions 
of X, 

N N 

-\7x{pf\7 + lY,yix-y^)>J2 ((1 - ^)«^(^ - y*) - f "'^(^ - y^^) ■ (50) 

i=l 1=1 

Proof. This follows immediately from the previous lemma, using translation invariance and the fact that 9r{x — 
yi) < 1 since all the balls are non-overlapping, by assumption. □ 

To apply this corollary, let l{p) be a smooth, radial, positive function with l{p) = for \p\ < 1, l{p) = 1 for \p\ > 2, 
and < l{p) < I in-between. For some s > let 

Xs{p) - l{sp). (51) 

Note that with this choice of Xsip) the corresponding h(x) = 1 — Xs{x) is a smooth function of rapid decay and hence, 
by simple scaling, the corresponding potential wji{x) satisfies (for R < const, s) 

R^ f R^ 

|'W'fl(a^)| < const. — and / \wr{x)\ d'^x < const. — (52) 

J s^ 

for some constants depending only on I. Moreover, if \yi — yj \ > 2R for all i ^ j, then 

^ 1 

'^wr{x - yi) < const. (53) 

2=1 

independent of x and N. Later we are going to choose R s qT^I'^ (cf. Eq. (|72(l '). 



B. A priori bounds 

For N\ ^2 — N, let ^'jv(-'f, Y) be a sequence of normalized wave functions, antisymmetric both in the X and in 
the Y variables. We assume that Ni/L^ gi and N2/L^ qi as — > oo, with q = Q\ + Q2- Let 71 and 72 denote 
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the reduced one-particle density matrices of 4* at {X, Y) for the X- and F-particles, respectively, with Tr 71 = Ni 
and Tr72 = N2. Moreover, let Pa/ denote the following spectral projection of the Laplacian with periodic boundary 
conditions on A, given by the integral kernel 

Pm{x,x')^^ ^ e^P<^-^'^ (54) 

pe(2Tr/i)Z3 
p|<(67r2M/i3)l/3 

for X, x' £ A. Note that, by scaling, Tr [Pm] does not depend on L, and 

lini ±Tt[Pm]^1. (55) 

In Lemmas Eland El we derive some bounds for sequences of wave functions satisfying certain energy bounds. These 
lemmas apply, in particular, to the true ground state - as shown in the previous section. We call these bounds a 
priori bounds. 

Lemma 5. Assume that, in the thermodynamic limit (N —> 00, L — *■ c» with Qi ~ Ni/L'^ fixed), there is a sequence 
of states n{X,Y) such that 

\imsup ^{^n\H\^n) < ^ (67r2)'/' ^ + g^/^j + Cag^ (56) 

L— too 5 L J 

for some C > independent of g. Then, for i = 1,2 (and with 71, 72 being the one-body density matrices), 

limsup -^Tr [7^(1 - Pn,)] < const, gy/ag^/^. (57) 



Proof. We immediately have the trivial lower bound for non- interacting fermions 



hminf -^^(M'jvl -Ax- Ay\^n) > ^ {Gn^f^'' \gl^' + gl^^] . (58) 



L3 

To prove Ineq. (|57|l . however, we need the following refinement of (|58|l . which is proved in Eq. (4.13)]: 

gl^^{l + const. Ci ) + gl^^il + const. Ci)] , (59) 



liminf ^(^I^jvl - Ax - AyI^^n) > - (Ott^)^^^ limsup 



where Ci = N.^ ^Tr [7^(1 — PjvJ] for i ~ 1,2. Using (|56|l as well as the fact that the interaction potential is assumed 
to be positive, we immediately obtain (|57|l . □ 

Our second a priori bound concerns the nearest neighbor distance among particles. For given points yi, . . . , yN-2 hi 
R'^, let Iii{yi, . . . , j/AT^) be the number of j/i's with the property that the distance to the nearest neighbor among the 
other j/i's is less than 2R. 

Lemma 6. Assume that there exists a C > 0, independent of g ~ N/L^ , such that 

^{^n\H\^n) <Cg^^^. (60) 

Then 

(«'iv|/fl(2/i,...,yAfJ|*w) < const. iV(i?3p)2/3. (61) 
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Proof. With 5i denoting the distance to the nearest neighbor, we have 

/fl(yi,...,yjv,)< (2i?)2^^. 
The result now follows from the operator inequality 

Y^-H < const. -A, (62) 

1=1 « i=i 

which holds on anti-symmetric wave functions of N2 variables, and is proved in [itI Thm. 5]. Note that we again 
used the fact that the interaction potential is positive, and hence the kinetic energy is bounded above by the total 
energy. □ 



C. Putting it together 

For a lower bound, we can neglect the interaction among particles of the same kind. That is, we use 

H > (^-Ax + I^XY^ + (^-Ay + I^XY^ . (63) 

We are going to bound both terms separately using the a priori bounds of the previous section. In the following, 
we are going to treat only the first term, the lower bound on the second term can be obtained in the same way by 
exchanging X and Y. 

First, we decompose the Laplacian into a high and a low momentum part, as follows: 

A = vr(p)v + V (1 - r(p)) V. 

For g = N/L^, let fcp — (Gtt^p)^/'^ be the Fermi momentum (for spinless fermions), and let 

r(p) = max h - , 



We claim that 



Ni „ ^5/3 

-V, (1 - r(ft)) V. > I {6n^f' (64) 

To show this, we use the argument in Let (f>i{x), i = 1, . . . , Ni^ be any set of orthonormal functions with support 
in the cube A of side length L. We want a lower bound to the expression 

/ \p\Hi~np))\Up)?dpp. (65) 

Note that ^,{p) = (27r)-3/2(e*p-i^|0.)^ ^ith ( • | •) denoting the inner product for functions on the cube A. Since the 
4>i{x) are orthonormal, we have that 



Y \Mp)\^ < {2TT)-^e'P-^\e'P-n = (27r)-3i:3 

Hence (|H5| is bounded below by the infimum of f p'^ (1 — T{p)) ^(p) (fip over all < S,{p) < {2'k)^'^L^ with J ^(p) d?p — 
Ni. Since |pp(l — r(j3)) is a monotone increasing function of \p\, the infimum is attained by ^(p) = 9{{Qt:'^ Ni/ L^Y^^ — 
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\p\), with 9 denoting the Heaviside step function. Now r{p) = for \p\ < {Gn'^Ni/L^Y^^ < (67r^gi)^/'^, and thus we 
arrive at H64(l . 

For the high-momentum part, we use that 

for any s < l/fcp, with Xs{p) defined in 1)51(1 . Hence, we can use Corollary ^ to get a lower bound on this term. In 
order to be able to apply this corollary, however, we have to make sure that the y^'s are separated at least a distance 
2R. Let Y C Y he the set of i/j's whose distance to the nearest neighbor is at least 2R. Note that, by definition, 
\Y\ — N2 ~ Ib{Y). We are going to neglect the interaction with j/j's that are not in the set Y, which can only lower 
the energy. Hence we obtain, for a given configuration of Y, 

J2 -v.r(p,)v, + - V'j) > (1 - ^^^1) E ^yi^i)^ 

i—1 i,j i—1 

with 

Wy{x)^ J2 {{l~e)aU{x~y,)-^WR{x-y,)), (66) 

depending on e, a, R and s. 

We are still free to choose U (x) . A convenient choice is 



Uix) 



3 - i?3) ^ for Ro < Ixl < i? 



otherwise, 

but any other choice such that |[/(a;)| < const. R^'^ for R^ Rq will do for our purpose. 

Now let 'i/ff{X, Y) be a normalized fermionic wave function. We can express the expectation value of Wyixi) 

as 



N 



'^Nj=J nyTr [7y VFy] dr, (67) 



where 



J \^N{X,Y)\^dX (68) 
and 7y denotes the one-particle density matrix of ^'Ar(X, y) for fixed Y, i.e., 

Jy{x,x') = — / 'i/N{x,X2, ■ ■ ■ ,XNi,Y)'ifN{x' ,X2, ■ ■ ■ ,XNi,yy d^X2 - ■ ■ d^XNi- (69) 

ny J 

Note that < 7y < 1 and Tr7y = iVi. Moreover, J ny dY = 1 and J nyjy dY = 71, the one-particle reduced 
density matrix for the X-particles. 

Let P be a projection operator, and let 7 denote any fermionic density matrix, which is an operator that satisfies 
< 7 < 1 and Tr7 — Ni. Let W± be two bounded positive semi-definite operators, and let W — W+ — For any 
(5 > 0, we have 

TT[-fW] = Tr [PW] + Tr [{j -1)PWP] 

+Tt [7 ((1 - P)WP + PW{1 - P) + (1 - P)W{1 - P))] 
> Tr [PW+] (1 - (5) - Tr [PW^_] (1 + 5) 

- (1 + S-') i\\W+\\ + \\W-\\)Tt [7(1 - P)] - ||VK!|Tr [P(l - 7)], 
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with II • II denoting operator norm. Now let P = Pni be the operator defined in 1)54(1 . and W = Wy- We choose W+ 
to be the terms in ((66|l containing U (x), and W- the ones containing wn{x). We then have, using J U{x) (fix — At:, 



T,[PW+]^^^I^ E U{x-y,)d^x 



N2 - Ir{Y) - const. — 



The last term in square brackets bounds the number of yj 's in Y that are at least a distance R away from the boundary 
of the box. Since the distance between the yj 's is bigger than 2R by assumption, the number of such yj 's close to the 
boundary is bounded by const. L^/i?^. By LemmaEl 

nylRiY) dY = {^N \Ir{Y)\ ^n) < const. N{R^gf/'' (70) 

if ^'jv(X, y) is an approximate ground state. As already noted in Eq. 1(55(1 . Tr [PjvJ can be replaced by A^i in the 
thermodynamic limit. 

Analogously, using 1(52(1 . we get an upper bound 

Tr[PW^_] < const.— 

Moreover, using l(55|) and the fact that the distance between j/j's contributing to Wy is at least 2R, we find that 

||W^yHoo< ||Ty+|| + ||W-|| < + const " 



The a priori bound in Lemma [S] implies that, for large enough iV, 

J nyTi [7y(l - P)] dY - Tr [71(1 - P)] < CNi^a'^Qf'^, (71) 

where 71 is the one-particle density matrix (for the X-particles) of any approximate ground state. The same bound 
is true for Tr [P{\ - 71)] = Tr [71(1 - P)] + Tr [P - 71], since iVf ^Tr [P - 71] ^ as iVi ^ cx3 (see Hence, 
collecting all the bounds, and applying the same arguments also to the second term in ((6H(I . we arrive at the lower 
bound 



lim 



^Snagig^ il-e-S- s^Gn^gf^^ - ] - Cag^iR^gf/"" 



es^R 



for some C > 0. 
We choose 



R = g-'/-' {ag'/'f^' , s = g-'/^ {ag'^') , £ = <5 = {ag'/') (72) 

and obtain, for small g, 



lim -^EoiN,,N2,L) > | (Gn^f' [pf ^ + gpl + Snagig^ - const, ag^/^^i/a^ Vis 
This finishes the proof of the lower bound. 
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VI. THE TWO-DIMENSIONAL GAS 

We now comment on the necessary changes in considering the 2D gas instead of the 3D gas. We start with the 
lower bound to the ground state energy. The analogue of Lemma 01 in 2D is the following lemma, which generalizes 
the corresponding result used for bosons in 2D in 9]. Its proof can again be found in the appendix. 

Lemma 7. For R > Rt^, let 6b{x) denote the characteristic Junction of a disc of radius R centered at the origin, i.e., 

Or{x) = 1 if \x\ < R and — otherwise. Let x{p) ^2 a radial function, < x{p) ^ 1; such that h{x) = (1 — x){^) 
bounded and integrable. Let 

Ir{x)^ sup \h{x - y) - h{x)\, (73) 
\v\<R 

and 

wr{x) = -fR{x) f fR{y)d\. (74) 
Let U{x) be any positive, radial function, supported in the annulus Rq < \x\ < R, with 

U{x)ln{\x\/a)d^x ^2t:. (75) 



Then, for any e > 0, 



B2 



'Vxip)ORix)x{p)y + ^v{x) > (1 - e)U{x) - J [(2^)- V [/(y) d^] WRix). (76) 



In the application, we choose, as in 



, iy(R)-^ for Ro < \x\ < R 
U{x) = <( _ I I _ ^^^^ 

otherwise, 



with v{R) determined by condition jTSJ, i.e., 

.R ^ 
v{R) = I \TL{r/a)rdr = — 
JRq 4 



T>2 p2 

RHn^-Rlln-^^ 
a^e a^e 



(78) 



Using that a < Rq < R we get the bounds 

i(i?2-i?o)2 (In i?/a-i) <i/(i?) < ^R^lnR/a, (79) 

from which, in turn, we get upper and lower bounds on J U (x) d^x = i'{R)^^it{R^ — Rq). 

Moreover, we again choose x{p) as in (|FT|) . with R < const, s. Inequalities (jE^ - lfS^ then have to be replaced in 
the 2D case by 

R^ f R^ 

\wr{x)\ < const. — and / (a;)| d^x < const. — (80) 



s4 



and 

N 



^ wr{x ~ yi) < const. (81) 

i=l ^ 

in case that \yi — yj \ > 2i? for all i ^ j. 
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The a priori bounds of Subsect . I V Bl can be obtained also in the 2D case. The proof of Lemma [S] works in the same 
way, with the appropriate changes in the expression of the kinetic and interaction energy, of course. For the proof of 
Lemma El we note that the analogue of the inequality (|62|l does not hold in 2 dimensions. However, a 'relativistic' 
version of it is true, namely that 

^ — < const. ^ v^-Ai (82) 

^=l °' i=l 

on antisymmetric functions of N2 variables j/i e . Ineq. H82f) can be proved in a similar way as the proof of H(j2f) in 
(irt. It implies that 

(*w|/i?(yi,...,y7vj|^'jv> <2i?Tr[V^7] < 2i?(TrhA7])i/2(Tr7)i/2 <const.iV(i?2^,)i/2 (33) 

in 2D, replacing H61|l . Here 7 denotes the one-particle density matrix (for the Y particles) of ^'jv(X, F), and we have 
used Schwarz's inequality as well as the assumption Tr [—A 7] < const. N g for an approximate ground state. 

With the a priori bounds in hand, we can proceed along the same lines as in Subsect. IV CI to obtain a lower bound 
to the ground state energy. The optimal choice of the free parameters e, 5, R and s in 2D turns out to be 

R=Q-'' \, , , s = ^ ^-A- \ 



|ln(a2£.)|3/20 ' |ln(a2£»)|i/20 ' | ln(a2e)|i/io " 

This yields the lower bound in Theorem [3 

Our last task is to derive the upper bound in Theorem |2J It turns out that obtaining this bound is actually much 
easier than in the 3D case. The reason for the rather complicated construction in 3D was the very small interaction 
energy ~ aq per particle, which forced us to choose the particle number in a box to be quite large, namely n ^ 1/(0"^ p), 
in order to have negligible finite size effects. This resulted in a trial wave function with very small norm. In 2D, 
however, it is possible to choose the particle number in each box much smaller, such that the norm of the trial wave 
function is close to one. If we take the analogous function as in l|ll|) - (|12(l . with s = 2i? and with Lp(x) now being the 
solution to the zero-energy scattering equation in 2D, cut off at an appropriate radius i?, then a simple bound as in 
the proof of Lemma|3| shows that 

k^Wl > 1 - const. n{R^g). (84) 

Hence we have to choose the box size £ and R such that nR^g^ ^ 1, with n ^ gi'^. Moreover, the restriction on 
having a negligible finite size effect is n^^^ '3> \ \n{a^g)\ (compare with l|19|) ). If we choose R = g~^^^\ \n{a^ g)]"" for 
large enough a, then all these conditions are easily fulfilled. In calculating the kinetic energy in H17|l and (|18|l . we can 
then just use the simple bounds g{x) < 1 and f{x) < 1. We demonstrate this on the analogue of the term II in H17|l 
in 2D. Namely, with ^(x) given as in (|25|l . 

J [\WxF{X,Y)\^ + ^vxYF{X,Yf] D,,{XfD,n{YfGn{XfG^{YfdXdY 

a^^-yJ)Dn(X)DUY)dXdY^ I gfjx)g^{y)ax~y)d\d\. (85) 

Here we have also used the fact that the integrand vanishes whenever two particles of the same kind are closer together 
then a distance s > 2R, in order for to hold. We can then proceed using Young's inequality on the last term, as 
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in H37() - (|39|l . The leading term from the interaction energy then comes from 

/ (|V^(x)p + i.(x)|^(.)n d'x = < -i^ fl + const, ^"'yf^,^' ) . (86) 

J\x\<R ln(i?/a) \ln{a^g)\ \ |ln(a^^?)| / 

The other terms in the upper bound can be treated in the same way. It turns out that, choosing a large enough, all 

the other error terms besides the one in H86|) are of lower order in g. We omit the details. This results in the upper 

bound stated in Theorem [3 



APPENDIX A: PROOF OF LEMM AS H AND □ 



We start with the three-dimensional case, Lemma 0] It suffices to show that the operator inequality H49|) holds for 
the expectation value with any smooth function ip{x) of compact support. Given such a ^{x), define the function 
^(x) by its Fourier transform ^(p) = x{p)i'{p)- We thus have to show that 

[|Ve(a;)|2 + i«(a;)|^(x)p]d3a;> / [(1 - £)aC/(a;)|^(x)|2 - rf'x. (Al) 

\x\<R. JE3 L ^ J 

Let ip{x) denote the solution to the zero-energy scattering equation |SJl, subject to the boundary condition 
lim|j.|^oo v(^) = 1- Let be a complex-valued function on the unit sphere §^ , with J^2 Wl'^ — 1- We use the 
same symbol for the function on taking values t^(a;/|a:;|). For 'ipi^) above, consider the expression 

A= iy{x)WC{x) ■^(p{x)(fx + ^ i v{x)'il){x)*Lp{x)v{x)d^x. 

J\x\<R J 

We note that the last integral makes sense even in the case when v{x) has a hard core; in this case, ^v{x)ip{x) has to 
be interpreted as the (non-negative) measure A(^(a;) (see Eq. Q). By using the Cauchy-Schwarz inequality, we can 
obtain the upper bound 

I^P < ( / [|V^(x)p + \v{xmx)\''] dPx] ( f [|V^(x)p + v{xf d^x] . 

\J\x\<R J \J\x\<R J 

Since ip{x) is a radial function, the angular integration in the last term can be performed by using Jg2 W\'^ = 1- The 
remaining expression is then bounded by a because of Jj^g (|V(p(x)p -I- ^w(x)|(p(x)p) d^x = Aira, as pointed out in the 
beginning of Section llVI Hence we arrive at 

[|VC(x)p + '^v{x)\^{x)\^] d\ > (A2) 

,x\<R. 

for any choice of v as above. It remains to derive a lower bound on \A\'^ . 

Note that (p{x) is a radial function with |V(y5(x)| = a/R^ for \x\ — R. Hence we obtain, by partial integration, 

v{x)\7C{x) ■'^'f{x)d\ ^ ~ ( C{x)'^{x)^f{x)d^x + ^ I C {x)v{x) dujR, 

x\<R J\x\<R ^ J\x\=R 

where dtOR denotes the surface measure of the ball of radius i?, and we used the fact that Vv{x) ■ Wip{x) = 0. Now, 
by definition of h{x), ^(x) = iIj{x) — {2'K)^^^^h * ^{x), where * denotes convolution, i.e., h * ipix) = J h{x — y)4'{y) d\j. 
Using the zero-energy scattering equation (j^J for 1^9(2;), we thus see that 

I r{x)v{x)dwR - (27r)-3/2 « f {h*,p)*(x)v{x)dwR 

^ J\x\ = R ^ J\x\=R. 

+ (27r)-3/2 f {h*',p)*{x)i^{x)A(p{x)d^x. (A3) 

J\x\<R 
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The last two terms on the right side of (|A3|) can be written as (note that /i is a real- valued function) 



(2^) 



-3/2 



h{y - x) dfi{y) 



(A4) 



where d/i is a (non-positive) measure supported in the ball of radius R. Explicitly, d^{y) = —aR '^i^{y)S{\y\ — R)d^y + 
iy{y)Aip{y)d^y. Note that J dii{y) — 0, and also J d\^{y)\ — 2a Jga 1*^1 < 2a\/An (by Schwarz's inequality). Hence 



h{y - x) dn{y) 



< 2aV^nfR{x), 



with fji{x) defined in 147|) . The expression HA4|I is thus bounded from below by 



in > -(27r)-3/22aV4^ J \ip{x)\fR{x) d^x > -a (^J \^p{x)\^WB.ix) d^x^ ' , 



(A5) 



where we used Schwarz's inequality as well as the definition of w_r(x) H48f) in the last step. Note that this last 
expression is independent of v(x). 

The only place where v{x) still enters is the first term on the right side of ljA3|l . By choosing v{x^ to be the 
restriction of "0(2;) to the sphere of radius i?, appropriately normalized, we obtain from (|A3p - (|A5l) 



A> 



R 



\ 1/2 

^(x:)\^dujR\ -ai I \ij{x)\''wR{x)d^x 



\x\ = R 



1/2 



Using again the Cauchy-Schwarz inequality, we see that, for any e > 0, 



\A\'>^{l-e) 



\x\=R 



\ip{x)\^ duoR - 



\'iP{x)\'^wr{x) d^x. 



(A6) 



In combination with (jA2|) this proves the desired result ljAl|l in the special case when U{x) is a radial (5- function 
sitting at a radius R, i.e., U{x) — R~'^S{\x\ — R). The case of a general potential U{x) follows simply by integrating 
this result (i.e., Ineq. (|A1I) for this special U{x)) against u{R)R^dR, with u{R) — U{x) for \x\ = R, noting that 
/ u{R)R^dR = 1 and that wii{x) is pointwise monotone increasing in R. 

The proof in the two-dimensional case. Lemma follows exactly the same lines. Note, however, that the solution 
to the zero-energy scattering equation can not be normalized by lim|3,|^oo ^i^) = 1 in 2D, but we can normalize it 
such that (p{x) — 1 for |a;| = R. It then follows that (p{x) — ln(|a;|/a)/ln(_R/a) for Rq < \x\ < R, and also that 
/|^l<^(|V^(x)p-f ii;(x)|cp(a;)|2)(i2a; ^ 27r/ln(i?/a) and J Aip{x) d'^x = 27r/ln(i?/a) (see the appendix in |g|). The rest 
of the proof is unchanged, with the result that 



[mx)\' + '^v{xM{x)\'] d'x>—^ 

x\<R m[n/a) 



(l-s) 



R 



\i;ix)\'du;H-- I mx)\'wR{x)d'x 

x\=R. £ 



(A7) 



instead of ljA2|l and (|A6|I . Multiplying this inequality with u{R)R\n{R/a), where u{R) = U{x) for |x| — R, and 
integrating over R using H75(l . we arrive at the desired result. 
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